Griffiths phase manifestation in disordered dielectrics 
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We predict the existence of Griffith phase in the dielectrics with concentrational crossover between 
dipole glass (electric analog of spin glass) and ferroelectricity. The peculiar representatives of above 
substances are KTaOs : Li, Nb, Na or relaxor ferroelectrics like Pbx-^La^Zr^^Tia.^Oz. 
Since this phase exists above ferroelectric phase transition temperature (but below that temperature 
for ordered substance), we call it "para-glass phase". We assert that the difference between para- 
electric and para-glass phase of above substances is the existence of clusters (inherent to "ordinary" 
Griffiths phase in Ising magnets) of correlated dipoles. We show that randomness play a decisive role 
in Griffiths (para-glass) phase formation; this phase does not exist in a mean field approximation. 
To investigate Griffiths phase properties, we calculate the density of Yang-Lee (YL) zeros of partition 
function and found that it has "tails" inherent to Griffith phase in the above temperature interval. 
We perform calculations on the base of our self-consistent equation for long-range order parameter 
in the external electric field. This equation has been derived in the random field theory framework. 
Latter theory automatically incorporates both short-range (due to indirect interaction via transversal 
optical phonons of the host lattice) and long-range (ordinary dipole-dipole interaction) interactions 
between impurity dipoles so that the problem of long-range interaction consideration does not appear 
in it. 



PACS numbers: 64.70.Pf, 77.80. Bh 

It was shown by Griffiths H that the free energy of a 
dilute Ising ferromagnet is a nonanalytic function of the 
external magnetic field for all temperatures between crit- 
ical value T c (x) (x is concentration of lattice sites without 
Ising spins) and T c (l) of the corresponding nondilute sys- 
tem. The manifestation of nonanalyticity in the original 
article of Griffiths was the distribution of the zeros of the 
partition function Z (the Yang-Lee (YL) zeros) in the 
plane of complex magnetic fields H — i6. Namely, be- 
low T c (l) YL zeros appear arbitrarily close to the point 
H = 0, implying a zero radius of convergence for expan- 
sions of thermodynamic quantities in powers of H. The 
" visible sign" of such nonanalyticity is the appearance of 
"tail" in density of YL zeros p{6) at T c {x) < T < T c (l). 

Dielectric systems are, in many respects, different from 
magnetic ones. The main difference is that main interac- 
tion in dielectrics is long-range dipole-dipole interaction, 
which is small relativistic correction to the short-range 
exchange interaction in magnets. It was shown in 
that former interaction in disordered dielectrics does not 
prevent the occurrence of aforementioned "tails" in p(9) 
and thus to realization of Griffiths phase analog. 

In the present paper we investigate the manifesta- 
tion of Griffiths phase in the dielectrics, where both 
long-range ferroelectric order and dipole glass (say, "di- 
electric spin glass") phase can occur depending on im- 
purity dipoles concentration. The peculiar representa- 
tives of such substances are KTaOs : Li, Nb, Na, 
where Li, Nb or Na, being off-center ions, constitute 
the impurity dipoles with some discrete number of per- 
missible orientations in host KTaOs matrix (see [p| for 
details). Other peculiar representatives of aforemen- 



tioned class of substances are relaxor ferroelectrics like 
Pb1-xLaxZro.Q5Tio.35O3 (PLZT), which can be consid- 
ered as some reference phase (in the case of PLZT it is 
PbZro.65Tio.35O3 having long-range ferroelectric order), 
"spoiled" by intrinsic defects of differnt kind, playing a 
role of impurity dipoles (see Q for details). 

Although the description of aforementioned disordered 
dielectrics is still not complete at present, we suggested 
so-called random local field method (see, e.g. j|] and ref- 
erences therein) for their description. This method cap- 
tures remarkably well the peculiarities of physical prop- 
erties of aforementioned substances, namely concentra- 
tional crossover between ferroelectric and dipole glass 
phases, existence of mixed ferro-glass phase and non- 
exponential long-time relaxation in glass and ferro-glass 
phases (|). 

In this method, we derive the equation for long-range 
order parameter L. This equation was derived self- 
consistently through the distribution function of random 
fields, acting between impurity dipoles (see [|| for de- 
tails). The specific form of this equation depends on 
particular number of orientations of impurity dipole in 
the host crystal matrix. For the simplest case of the im- 
purity dipole with only two permissible orientations this 
equation has the form 

tanh (f3E) f(E, L)dE, 
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f(E,L) = ±- [ exp(iEp) x 
x exp [Fi(p] + iF 2 (p)j dp, 
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where n is impurities concentration, f(E, L) is a distri- 
bution function of random fields (see, e.g. ||), Fi(p) is 
"responsible" for dispersion of random fields, Fz{p) - for 
mean value of random field, K, zz is a component of tensor 
of interaction JC a f3(r) (a,/? = x,y,z, f is interimpurity 
separation) between impurity dipoles in highly polariz- 
able dielectric host . General form of K. a p reads (see, e.g. 
H and references therein) 
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where eo is host lattice static dielectric permittivity, 
d* is the effective dipole moment of impurity (see, e.g. 
||), r c is host lattice correlation radius (two impuri- 
ties at a distance less or equal r c "feel" each other). 
The physical reason for appearance of r c and functions 
fi(x), f2(x) oc x 4 exp(— x) is indirect interaction of impu- 
rity dipoles via host lattice soft phonon mode (see, e.g., 
H and references therein). If latter interaction is ab- 
sent, we have r c — * 0, fi(r/r c ), /2(r/Y c ) — > so that (||) 
gives ordinary dipole-dipole interaction. It is seen that 
interaction (|J) incorporates both long-range part (term 
(3m a m/3 — 5 a f3)/r 3 ) and short-range part (terms contain- 
ing fi(x) and f2(x)). The explicit form of the functions 
fi and /2 can be found, e.g. in ^,^|. 

It can be shown (see, e.g. Q), that at high impurity 
concentrations distribution function of random fields has 
Gaussian form. This form can be obtained from (|l|) by 
the expansion of integrands in Fi^(p) U P to first leading 
order. This gives 
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Since the case of Gaussian distribution function implies 
the strongest (both long-range and short-range) interac- 
tion between impurity dipoles (and consequently their 
clusters) this case can be regarded as a most difficult 
for Griffiths phase to occur. We shall demonstrate the 
manifestation of Griffith phase for this case with under- 
standing that for lower impurity dipoles concentration 
the realization of Griffiths phase is easier. 

The equation for order parameter for Gaussian distri- 
bution function assumes the form 
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x sin(pE) sin p(£ + E L)dpdE. 
Introduce following dimcnsionlcss variables 

£ E kT 

— = h, pE = y, — = x, —— = t, 



(4) 



f 



E 2 



I57TZ 



(5) 



In these variables the equation for L assumes the form 
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x sin (xy) sin y (h + L) dxdy. 
Integration over x in gives 

y 2 \ siny(h + L) 
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We shall investigate the Griffiths phase properties on the 
base of Eq (0). 

The equation for phase transition temperature r c can 
be obtained from ([?]) at L — > (of course, at h = 0) 

'-I ' 7^—) . , t I- (8) 
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The equation for critical concentration of impurities can 
be obtained from (||) at r c — > 0. It reads 
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Equations (||J9|) determine the phase diagram of the sys- 
tem under consideration. It is shown on the Fig. I. 
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FIG. 1. Phase diagram of the substances under considera- 
tion. 
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The theorem of Yang and Lee states that zeros of a 
partition function Z lie on the unit circle h = i9 in the 
plane of complex variable z = exp(—2(3H), (/3 = l/k^T) 
JjJ. To calculate the density of YL zeros, we put in (Q) 



h — i9 1 L = L% 
We have from (^) 
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The equations (0),(|l^) are the main theoretical result 
of this paper. To find p(9) from ([[§), we should find (at 
given temperature and concentration) L\ and L 2 from 
(pi]). GrifHts phase is realized in the substances under 
consideration if p(9) has "tails" near 6 = at r c < r < 1. 
This is because in the " language" of equations (||) , @ the 
phase transition temperature of nondilute system T c (l) 
corresponds to the quantity Eq, determining the phase 
transition temperature in a mean field approximation 
T C MF, i-e. T c (l) = E = T c mf- 
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inherent to Griffiths phase for dilute system. For nondi- 
lute (completely ordered system) there is no Griffiths 
singularities. This clearly demonstrates the occurrence 
of Griffiths phase in the dielectrics with concentrational 
crossover between dipole glass and ferroelectricity. Fig. 2 
permits us to show the region of existence of Griffiths 
phase in the phase diagram of substances under consid- 
eration. 



Since Griffiths phase occurs at temperatures higher 
than that of ferroelectric phase transition of disordered 
system, but lower than that of the ordered (nondilute) 
system, we can call this phase "para-glass". Glassy 
behaviour is due to clusters of impurity dipoles (see 
p|,p|Jl^J^|). The manifestation of these clusters is indeed 
"tails" in p(0). It is seen from Fig. 2 that at temperatures 
than T c mf there is no more "tails" in p(9) and para-glass 
phase becomes conventional paraelectric phase. 

Other peculiar feature of all curves is the existence 
of some threshold value 9 e , below which there are no 
YL zeros (p(9) = 0). This is so-called YL edge, which 
is manifestation of ordered state of the system (para- 
magnetic phase of nondilute system, see PJnffl ) . The 
"tail-like" approach to 9 e at smaller z is due to disorder 
(glassy effects) in the system. Thus at sufficiently small 
z (z = 0.05 in Fig. 2) the system is in a mixed ferro-glass 
phase, exhibiting both "order" and "disorder" features. 
The difference between ferro-glass and para-glass phase 
is that while in former the overall long-range order exists 
in the system, in latter it does not exist. Our supposi- 
tion that para-glass phase corresponds to superparamag- 
netism, well-known in spin glasses. 



To find the equation for 9 e we put L\ = in flll] 
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FIG. 2. Density of YL zeros. Characteristic "tails" inher- 
ent to Griffith phase are shown. Curve labeled "MF" is also 
realized in mean field approximation (see Eq. (|Ts|) ) . 

We plot the dependence p(9) for different r and z in 
the Fig. 2. First of all, it is seen the presence of "tails", 



The dependence e (r) is shown in Fig. 3. It is seen that at 
z > z cr = 1/15 (ferroelectric and ferro-glass phase) 9 e = 
at t = t c , while at z < z cr 9 e ^ for all temperatures. 
To trace the transition from ordered to Griffith phase it is 
instructive to investigate the concentrational dependence 
of 9 e at zeroth temperature. 
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FIG. 3. Temperature dependence of 8 e at different impu- 
rity concentration. Curve labeled "MF" is also realized in 
mean field approximation (see Eq.(po[)). 
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Integral in (|15|) can be calculated analytically giving 
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Denoting /i = 9 + L 2 , we obtain 
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The dependence ( |17| ) is reported in Fig. 4. It can be 
shown that equation 9 e (X) = is satisfied at z — z cr = 
1/15. 



FIG. 4. Concentrational dependence of 
perature. 
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Let us finally show that para-glass phase cannot be re- 
alized in the ordered dielectric. To do so we will demon- 
strate the absence of " Griffiths tails" for nondilute sys- 
tem. This corresponds to mean field limit (z — > oo) in 
random field method. In this limiting case we have from 
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This equation is indeed the mean field equation for or- 
der parameter in Ising model. In this case t c mf = 1- 
Substitution © to <M) gives 
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To calculate 9 e (r) in this case we should put L\mf = 
in (flj). This after some algebra yields 
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The dependences p(9) and 9 e (r) for nondilute system are 
shown in Fig. 2 and 3 (curves, labeled "MF") respectively. 
It is seen the absence of "tails" in p{9). This shows, that 
random fields (via their distribution function f(E, L)) 
play major role in Griffiths phase formation. 

In the present paper we have shown the existence of 
Griffiths phase, inherent to disordered Ising magnets, 
in the dielectrics with concentrational crossover between 
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dipole glass and ferroelectricity. We call this phase " para- 
glass" because it combines the paraelectric phase features 
(existence of YL edge 9 e ) and "disorder features" (tails 
in p(9)). This phase resembles very much superparam- 
agnetic phase, peculiar for many spin glasses. 

Since physical properties in glass and ferro-glass phases 
of aforementioned substances are different from those in 
para-glass (Griffiths) phase, we can expect some pecu- 
liarities of latter phase, which might be helpful in its 
experimental observation. One of such features may be 
nonexponcntial relaxation of polarization, different from 
that in glass and ferro-glass phases. Its investigation is 
interesting problem and can be done within above for- 
malism. 

Unfortunately, we don't know the direct experimental 
observation of the Griffiths phase in the aforementioned 
substances. For the description of possible (direct or in- 
direct) experiments our approach can be easily extended 
to the case of nongaussian distribution function of ran- 
dom fields as well as for tunneling of the aforementioned 
off-center impurities between their permissible orienta- 
tions (i.e. between minima of their multi-well potentials) 
in the host dielectric. 
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